Quantitative study of the violation of /^-factorization 
in hadroproduction of quarks at collider energies 



Hirotsugu Fujii a , Frangois Gelis b and Raju Venugopalan c ' d 
a Institute of Physics, University of Tokyo, 
Komaba, Tokyo 153-8902, Japan, 
b CEA/DSM/SPhT, 91191 Gif-sur-Yvette cedex, France, 
c Physics Department, Brookhaven National Laboratory, 
Upton, N.Y. 11973, U.S.A. 
d ECT-k, Villa Tambosi, Strada delle Tabarelle 286, 
1-38050, Villazzano(TN), Italy. 

(Dated: February 2, 2008) 

We demonstrate the violation of fcx-factorization for quark production in high energy hadronic 
collisions. This violation is quantified in the Color Glass Condensate framework and studied as a 
function of the quark mass, the quark transverse momentum, and the saturation scale Q s , which 
is a measure of large parton densities. At x values where parton densities are large but leading 
twist shadowing effects are still small, violations of fcx-factorization can be significant - especially 
for lighter quarks. At very small x, where leading twist shadowing is large, we show that violations 
of fcx-factorization are relatively weaker. 



The k± -factorization formalism devised originally 
for heavy quark production in hadronic collisions, sys- 
tematically resums powers of a s ln(s/g^) in perturbative 
QCD. These contributions are important at high energies 
when the transverse momentum q± of the hnal state is 
much smaller than the center of mass energy y/s. In the 
k± -factorization framework, the quark production cross- 
section is expressed as the convolution of a hard matrix 
element and distribution functions from each of the two 
hadrons. These "unintegrated" gluon distributions de- 
pend on the longitudinal momentum fraction x and the 
transverse momentum k^ of the gluon taking part in the 
hard scattering. When integrated over all transverse mo- 
menta up to some hard scale M 2 , these distributions give 
the more familiar gluon distribution x G(x, M 2 ). 

Even though logarithms due to gluon branchings are 
resummed in the framework of refs. only one hard 
gluon from each projectile participates in the reaction. 
It is interesting to consider whether the fcx-factorization 
framework can be extended beyond the single hard scat- 
tering (leading twist) case. The simplest experiments for 
studying hard multiple scattering (higher twist) effects 
are those where one of the projectiles is dilute and the 
other is dense. These include a) p-p collisions in the for- 
ward/backward fragmentation region at the Relativistic 
Heavy Ion Collider (RHIC) and the Large Hadron Col- 
lider (LHC), where large x's in one projectile and small 
x's in the other are probed, and b) proton or Deuteron 
collisions off large nuclei at RHIC and LHC. Understand- 
ing the validity of k± factorization in this multiple scat- 
tering regime is important for a quantitative understand- 
ing of final states at RHIC and LHC. 

The Color Glass Condensate (CGC) 2], wherein the 
two projectiles are sources of classical color fields, is a 
powerful framework to study multiple scattering effects. 



The leading twist results of refs. [lj for quarkpair pro- 
duction are easily recovered in this formalism 3] keeping 
only terms that are of lowest order in the sources, fcx- 
factorization was studied in the forward p-p/p-A case, 
where one keeps the lowest order in one of the sources 
and all orders in the other source - this power count- 
ing is discussed further after eq. I@J. In refs. 0,0) this 
formalism was used to obtain the cross-section for gluon 
production in p-A collisions, which was found to be fcj_- 
factorizable in 5, 6]. Quark production in p-A collisions 
was computed in the CGC formalism in 0, 0- R was 
shown explicitly in that both quark pair produc- 
tion and sin gle i nclusive quark cross-sections are not fcj_- 
factorizable [l2|. In this note, we quantify, for the first 
time, the magnitude of the violation of fcj_ factorization. 

The generalization of the leading twist formula of 
refs. Q for the inclusive quark production cross-section 
in p-A collisions gives pfl, 
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where the shorthand notation corresponds to the explicit 
expressions, 
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is the well known Lipatov effective vertex. In eq. , 
feij^ and fc2_L are respectively the transverse momenta 
transferred from the proton and the nucleus, and p_|_ 
must be understood as + k 2 ± — q±. All the relevant 
information about the proton and the nucleus is encoded 
in the function ip v and in the various <fi A 's respectively. 
These are defined as 
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for the nucleus. R p and R A are the radii of the pro- 
ton and the nucleus respectively. U is a Wilson line in 
the adjoint representation, and U is a Wilson line in the 
fundamental representation. The function (p p , expressed 
here in terms of the color charge density pi in the proton, 
is the usual leading twist unintegrated gluon distribution 
of the proton |3j. The functions <f> A , on the other hand, 
include higher twist re-scattering corrections to all or- 
ders in /?2 Aj_> the ratio of the color charge density in the 
nucleus divided by the square of the momentum l\_ trans- 
ferred from the nucleus. These correlators, and models 
to compute these, will be discussed further shortly. We 
will only mention here that the three 4> A 's obey, in full 
generality, the sum rules, 
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Equation has another important property. When 
integrated over p + , the Dirac traces that appear in this 
formula become independent of q + , or equivalently, of 
the rapidity of the produced quark. In a model like the 
McLerran-Venugopalan (MV) model ^(|; where the dis- 
tributions ip p and (j) A are independent of x, the single- 
quark cross-section is rigorously rapidity independent. 
When the distributions have an x dependence, it is a pri- 
ori not legitimate to integrate the Dirac traces over p + 
because of the p + dependence in the x variables. How- 
ever, doing so is a good approximation if the p + depen- 
dence of the Dirac traces is much stronger than that of 
the unintegrated gluon distributions. In this approxima- 
tion, the integral over p + of the Dirac traces can then be 



done in closed form. One of these is particularly simple 
and gives 
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Therefore, the sum rule obeyed by the distribution (j> q ^ q 
applies for this term and we can replace the factor 
(C p /N)<j% q by § q f without changing the result. Thus 
the only term in eq. Q that truly breaks fcj_-factorization 
is the term involving the 3-point correlator cj) qq ' 9 . For the 
other two terms in the cross-section, the result of the an- 
alytic integration over p + is rather involved and will not 
be quoted here. 

We shall now compute the correlators in eqs. J3J and 
IpTjl in models which give good results in specific kine- 
matic regions. As discussed more extensively in , stud- 
ies of Gold-Gold and Deuteron-Gold collisions at RHIC 
suggest that the classical MV model, where the expec- 
tation values of the Wilson line operators are computed 
with Gaussian source functionals, is a good model for 
source distributions of large nuclei at moderately small 
x (x ~ 10 -2 ). In this kinematic region, where small x 
quantum evolution is not dominant, multiple scattering 
contributions can be interpreted entirely as higher twist 
effects. These die away rapidly with transverse momen- 
tum. When small x quantum evolution becomes signifi- 
cant, the Gaussian distribution of sources breaks down. 
The evolution of the source distributions at smaller values 
of x is described by a functional renormalization group 
(RG) equation Q which remains to be solved. However, 
at very high parton densities, color charges are strongly 
screened allowing the application of mean field techniques 
analogous to the well known Random Phase Approxima- 
tion. In this limit, the source functional is remark ably 
also a Gaussian, though (crucially) a non-local one 
Unlike the local MV model, which does not include non- 
trivial screening and recombination effects, one obtains 
leading twist shadowing from this Gaussian mean field 
model (GMF). As discussed in 0, it is conceivable that 
this regime might already be reached for x ~ 1(V 3 in 
large nuclei. In the absence of a full solution of the RG 
equations, the MV model and the GMF model corre- 
spond to two extremes of high parton density behavior. 
We expect more quantitative results for single quark pro- 
duction to lie between these two extremes. 

In Appendix A of ref. Q , we derived analytic expres- 
sions for the correlators of Wilson lines that appear in 
the <p A 's when the distribution of color sources in the 
nucleus is a Gaussian, be they local or non-local. We 
shall use these results here. The 2-point correlators for 
these models in momentum space are easily computed 
since they reduce to a 1-dimensional integral. This fol- 
lows from rotational symmetry in the transverse plane. 
The 3-point function 4> qq ' 9 is more computationally inten- 
sive because it is determined by a 3-dimensional Fourier 



3 



integral. fcx-factorization requires that the 3-point and 

2- point correlators are related by [7j 

<f>f'°(l ± ; k ± ) = i(27r) 2 [S(k ± ) + 5(k ± - l ± )] r/(l±) ■ 

(7) 

This can be checked by substituting this relation in 
eq. JIJ. The ratio, for the MV model, between the exact 

3- point function and the 2-point function <f> 9 ^ 9 as a func- 
tion of k± at fixed l± is displayed in figure ^ At large 
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FIG. 1: The ratio <j>f ' 9 {l ± - k ± )/^ 9 {l ± ) as a function of k ± 
(along a line parallel to l± in the k± plane) for various values 
of fx and Q\ = 2 GeV 2 . The thin lines represent the same 
ratio evaluated with the large N approximation for the 3-point 
correlator in the numerator. 

enough l±, this ratio indeed has two peaks centered at 
k± = and k± = lj_ respectively. The width of these 
peaks is roughly of the order of the saturation momentum 
in the nucleus, Q s . Eventually, when l±_ decreases below 
a value of order Q s , the two peaks merge into a single 
maximum centered at l±/2. This behavior of the three 
point correlator suggests that fcx-factorization is a good 
approximation if the produced quark is characterized by 
some large momentum scale (m or q±). Indeed, if this is 
the case, the typical k2± in eq. (JTJ is large compared to 
Q s , and eq. (JJJ is a good approximation of the two-peak 
structure of the exact 3-point function. 

In the large N limit, the 3-point function in the MV 
and GMF models is a product of two 2-point functions. 
The validity of the large N computation of the 3-point 
function (represented by thin lines in figure ^ is rea- 
sonably good. In figure |3 we compare the quark cross- 
sections obtained in the full calculation with those in the 
large N approximation computed for N — 3. They are 
close to one another since terms neglected in the large 
TV limit are of order 1/N 2 . We shall therefore perform 
all subsequent numerical calculations in the large N ap- 
proximation. 

We shall now discuss the violation of fcx-factorization 
in cross-sections. In Fig. |3 we plot the ratio of the com- 
plete result to the k± factorized result (obtained by using 
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FIG. 2: Comparison of the cross sections obtained with the 
exact 3-point function and with the large N approximation, 
for Q 2 a =2 GeV 2 , m = 1.5 GeV and m = 4.5 GeV. 

eq. Q for the 3-point function) as a function of the trans- 
verse momentum in the MV model, respectively for quark 
masses m — 0.15 GeV, m = 1.5 GeV and m = 4.5 GeV . 
The various curves in each figure correspond to different 
values of Q s . Extrapolating from saturation models of 
the HERA data (see for a discussion), one estimates 
Ql = A 1 / 3 (x /x) x , where x = 3 • 1(T 4 and A w 0.3. 
The values of Q 2 S shown in Fig. [3] correspond to a wide 
range of x. The values of Q\ = 1 GeV 2 and 4 GeV 2 
correspond to the central regions at RHIC and LHC re- 
spectively. Larger values of 15 and 25 GeV 2 correspond 
to very forward kinematics-between 4 and 6 units of ra- 
pidity in the proton direction in future proton-Lead ex- 
periments at the LHC; large values of Q s may also be 
accessed in proposed upgrades at RHIC. 

From the figure|3 one can clearly deduce the following: 

(1) The breaking of fcx-factorization is quite sensitive 
to the quark mass. The magnitude of the breaking is 
systematically smaller for heavier quarks because they 
are less sensitive to re-scattering effects than light quarks. 

(2) The magnitude of the breaking for the heavier 
quarks becomes significant only at values of the satura- 
tion scale that correspond to forward rapidities. Isolating 
these effects from other uncertainities in these measure- 
ments is difficult but can be done in principle. One might 
have to look simultaneously at other final states, such 
as hard diffractive final states, that are also sensitive to 
three point functions. 

(3) The magnitude of its breaking is maximal for 
q± ~ Qs- One recovers fcx-factorization when the quark 
transverse momentum is much larger than Q s . 

(4) If Q s remains smaller than the quark transverse 
mass, the breaking of fcx-factorization enhances the 
cross-section because re-scatterings push a few more QQ 
pairs above the pair production threshold. Conversely, if 
Q s is large relative to the quark mass, the cross-section is 
reduced at small transverse momentum because multiple 
scatterings typically transfer large transverse momenta 
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to the quark. 
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FIG. 3: Breaking of fcx-factorization in the MV model, re- 
spectively for strange, charm and bottom quark production. 

As discussed previously, these conclusions in MV 
model are valid when parton densities are high but quan- 
tum evolution is not significant. Computing the same 
quantities in the GMF model provides a good idea of the 
direction and magnitude of the effects when small x quan- 
tum evolution becomes significant. The results in figure 
0] for the production of charm quarks, display the same 
trends outlined previously. The magnitude of the break- 
ing of fc_L-factorization is relatively smaller in this model. 
Similar results were obtained in our previous study |5j 
of re-scattering effects in gluon production. An intuitive 
way to understand this result is to note that recombina- 
tion/screening effects lead to a depletion in the number 



of gluons available to re-scatter. For comparison, we also 
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FIG. 4: Breaking of fcx-factorization for charm quarks in the 
non-local Gaussian effective model of |Tll| . The dashed line 
shows the result (for Q 2 S = 15 GeV 2 ) in the MV model for 
comparison. 

plot the result for charm quarks at Q 2 S — 15 GeV 2 in the 
MV model. We expect results of detailed RG computa- 
tions to lie within the band outlined in Fig. 0] 
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